Schwinger boson theory of anisotropic ferromagnetic ultrathin films 
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Ferromagnetic thin films with magnetic single-ion anisotropies are studied within the framework 
of Schwinger bosonization of a quantum Heisenberg model. Two alternative bosonizations are 
discussed. We show that qualitatively correct results are obtained even at the mean-field level of the 
theory, similar to Schwinger boson results for other magnetic systems. In particular, the Mermin- 
Wagner theorem is satisfied: a spontaneous magnetization at finite temperatures is not found if 
the ground state of the anisotropic system exhibits a continuous degeneracy. We calculate the 
magnetization and effective anisotropies as functions of exchange interaction, magnetic anisotropies, 
external magnetic field, and temperature for arbitrary values of the spin quantum number. Magnetic 
reorientation transitions and effective anisotropies are discussed. The results obtained by Schwinger 
boson mean-field theory are compared with the many-body Green's-function technique. 



PACS numbers; 75.10.,]m, 75.30.Gw, 75.70.Ak 
I. INTRODUCTION 

Schwinger boson theories are very successful in-describ- 
ing magnetism in various quantum systems .ETlI Quah- 
tatively correct resuhs are mostly obtained even in the 
mean- field approximation, e.g.^ quantum Heisenberg fer- 
romagnets and antiferromagnets are well described for 
arbitrary spin quantum number for two and more spa- 
tial dimensions. In particular in two-dimensional (2D) 
isotropic Heisenberg magnets no spontaneoua-,njiagneti- 
zation is found for non-zero temperatures^uo in ac- 
cordance with the Mermin- Wagner theorem.B There are 
basically two reasons why Schwinger boson mean-field 
(SBMF) theory works well even in low dimensions: (a) 
Since the bosonized spin degrees of freedom are inte- 
grated over in the functional integral, spin fluctuations 
are taken into account and (b) the approach does not 
constitute an expansion around an^-ordered state, un- 
like Holstein-Primakov bosonization,H and thus works for 
both ordered and disordered ground states. 

The first aim of the present paper is to test the appli- 
cability of SBMF theory to another class of quantum- 
magnetic models, namely the Heisenberg ferromagnet 
with single-ion (on-site) anisotropies. We will restrict 
ourselves to the two-dimensional case, which is the most 
interesting one, but the same methods can be used in 
any number of dimensions. The anisotropic Heisenberg 
model in two dimensions is of particular interest, since 
it describes ultra-thin ferromagnetic films, including the 
important cases of the 3d transition metals iron, cobalt, 
and nickel. Thin ferromagnetic films are of great techno- 
logical relevance, e.g., for magnetic disk drives. They also 
show many interesting physical effects, of which we only 
mention the observed reorientation transitions-as-a func- 
tion of both film thickness and temperature. E2rE3 Much 
theoretical work has already been done on anisotropic 
films, which we briefly review below. However, it would 
be desirable to have an alternative approach for compar- 
ison. The second aim of this paper is to provide such an 



approach using Schwinger bosonization. We apply this 
method to the calculation of the magnetization with and 
without an external magnetic field. In addition we deter- 
mine effective anisotropies, which are crucial for compar- 
isons of measurements performed at finite temperatures 
T to ab-initio calculations of anisotropies at T = 0. The 
effective anisotropies describe the dependence of the free 
energy on the polar angles {9, (j)) of the magnetization. 

From the theoretical point of view, two-dimensional 
anisotropic ferromagnets are of particular interest, since 
the anisotropies reduce the symmetry of the system. For 
example, an easy axis leads to a discrete, Ising-type sym- 
metry, while a easy plane leads to an XK- type, symmetry. 
According to the Mermin- Wagner theoremja a sponta- 
neous magnetization at finite tempemtures is possible in 
the first case but not in the second.EJ We will show that 
SBMF theory yields qualitatively correct results for both 
of these cases. 

We now summarize previous work on thin magnetic 
films. So far, most approaches rely on simple mean-field 
approximations, which are of limited validity for low- 
dimensional systems. In particular, they do not satisfy 
the Mermin- Wagner theorem. A finitepiaagnetization can 
be induced by magnetic anisotropies, Ell which, however, 
cannot be described by a perturbative expansion around 
the symmetric state. Thus, improved approaches are 
needed, which should also allow to determine the mag- 
netic properties in the whole range of temperatures and 
for systems with nonequivalent lattice sites. For exam- 
ple, magnetic reorientation transitions have been stud- 
ied within spin-wave thjeories. In particulp:. Holstein- 
Primakov bosonizationO has been applied,E3 which is, 
however, limited to low temperatures. 

The magnetic properties of 2D anisotropic ferromag- 
nets have also been studied using many-body Green's 
functions, which allow cajlculations in the whole temper- 
ature range of interest .cJ^Zl This metliOjciJ|ias also been 
applied to thin films of several layers.EfTEj Within this 
theory the higher-order Green's functions are approxi- 
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mated by the so-called Tgablicov (or random-phase ap- 
proximation) decoupling^ i.e., one of the spin oper- 
ator components is replaced by its expectation value, 
whereas the expectation values of the other two com- 
ponents vanishes. Note that within this approach the 
spin commutation relations are preserved. By compari- 
son with a recent Quantum-Monte-Carlo calculation for 
an isotropic 2D magnet with S =i^-i/2 on a square lat- 
tice in an external magnetic field,ol2l it has been shown 
that this approaclu-yields a reasonable description for 
the magnetization.Ej In addition, the long-range mag- 
netic dipolc interaction can be treated as well within the 
same decoupling scheme. On the other hand, single-ion 
anisotropies have to be treated differently, since the iit. 
elusion of such terms leads to an ambigaus description.E3 
Anderson and Callen,LZl as well as LinesE3 have proposed 
interpolation schemes, which are valid only for small 
single-ion interactions as compared to the exchange cou- 
pling. For most cases of interest this condition is fulfilled. 
In order to calculate the effective anisotropies K,i{T) and 
the temperature-dependent magnetic reorientation, non- 
vanishing expectation values of at least two spin operjap. 
tor components have to be determined simultaneously.EZI 
Whereas the inclusion of the magnetic dipole coupling is 
straightforward, the inclusion of single-ion anisotropies 
turns out to be tedious, in particular for spin quantum 
numbers S' > 1, and for higher-order anisotropies. We 
emphasize that within these Green's-function approaches 
the magnetic anisotropies cannot be handled as small 
perturbations. Also, this approach suffers from the am- 
bigous decoupling schemes for these single-ion terms. 
Thus, a theory not suffering from these problems would 
be desirable. 

Finally, we summarize other pappcoaches. First, mean- 
field theories have been usedJ^^EJ which are particu- 
larly simple to apply for complicated systems, e.g., with 
nonequivalent atomic layers in a thin film or several non- 
equivalent lattice sites. Classical as well as quantum- 
mechanical Heisenberg spins have been studied allowing 
for an arbitrary direction of the magnetization. The 
inclusion of single-ion anisotropies is quite simple. A 
perturbative treatment of the anisotropies yields- |ajia^ 
lytical expressions for the effective anisotropies .CjO'l2I 
However, these approaches completely neglect both ther- 
mal and quantum fluctuations, which-are crucial for low- 
dimensional systems. Kawazoe et alE^ have treated the 
magnetic reorientation of thin ferromagnetic films within 
a continuum approach for the system and mean-field- 
like expressions for the magnetization density. A con- 
tinuum approach to the magnetic reorientation has also 
been used within a rcnormalization-group trcaJ|iacnt.cZI 

On-4;he other hand, following earlier worlO'Ej Millev 
et a/Jlj treat the magnetization of a thin film as a macro- 
scopic, classical magnetic moment and discuss its ground 
state as a function of anisotropy and external magnetic 
field. The anisotropy constants appearing in their work 
are to be understood as effective anisotropies obtained 
from experiment or from some more advanced calcula- 



tion. This approach has to assume that the effects of 
thermal and quantum fluctuations, the lattice structure, 
and the spin quantum number S can be described by 
a suitable renormalization of the anisotropy constants. 
From a theoretical point of view, the problem has just 
been transferred to the calculation of these effective 
anisotropies. This has to be done using more advanced 
methods such as Schwinger bosonization (discussed be- 
low) or many-body Green's functions. 

Finally, Monte-Carlo simulations have been performed 
for the effective anisotropies and magnetic reorienta- 
tion of ferromagnetic monolayers. So far, only classical 
spins have been studied for the anisotropic case,E3 while 
Quantjun-Monte-Carlo results exist for the isotropic case 
only.Ejfl More complicated systems such as thin films 
with several layers or including the dipole interaction 
have not been studied yet. 

Following the goals set above, we derive in Sec. || 
SBMF theories for films with second- and fourth-order 
single-ion anisotropy. We employ two different bosoniza- 
tion schemes, known as SU(A^) and 0{N) bosonization. 
In Sec. [II we show representative results and compare 



them with results from simple mean-field theories and 
from many-body Green's-function met hods. We focus 
on two areas of current interest: In Sec. HI A we discuss 
the change of magnetization as a function of external 
field and of temperature, in particular the reorientation 
tran sition s, and present a number of phase diagrams. In 
Sec. [II B we calculate temperature-dependent effective 
anisotropies. A summary and conclusions are given in 
Sec. ITvl. 



II. SCHWINGER BOSON THEORY 

In this section we summarize the SBMF theory for an 
anisotropic 2D ferromagnet, using two different, but re- 
lated, bosonizations. i-iThe first one employs the SU(2) 
symmetry of the spinas in the absence of anisotropies and 
leads to an SU(7V)-symmetric mean- field theory, where 
N oo. The second makes use of the mapping be- 
tween the Lie groups SU(2) and 0(3) and results in an 
0{N) mean- field theory.Q We first sketch the SU(iV) the- 
ory, which is described in greater detail in App. The 
derivation for the 0{N) theory is similar and we only 
present the results. 

The main idea of SU(2) Schwinger boson theory is to 
map the spiiirOperators onto boson operators 6| and 6| 
according toou 
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This mapping is exact; in particular the spin commu- 
tation relations are correctly reproduced. However, the 
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bosonization introduces unphysical degrees of freedom, 
which have to be removed by the constraint 



b\bT+blbi = 2S, 



(4) 



where S is the spin quantum number. Whereas the above 
scheme employs the SU(2) spin symmetry group, one 
can also take advantage of the mapping between SU(2) 
and the group 0(3) — both have the same algebra up 
to the choice of representation — to formulate an alter- 
native 0(3) Schwinger boson theory.El For both groups 
well-defined expansions around a mean-field theory are 
obtained by generalization to N boson fields and expand- 
ing in small 1/Ns In particular, in the limit N ^ oo 
mean-field theory becomes exact. 

Read and SachdevO have derived SU(iV) and 0{N) 
SBMF expressions for the magnetization of isotropic 2D 
ferromagnets, where their main interest was in the-quan- 
tum Hall ferromagnet at filling factor v = l.tilH Re- 
cently, the leading-order fluctuation corrections, i.e., the. 
first order of the 1/iV expansion, have been calculated ^ 
and compared to Quantum-Monte-Carlo simulationaSJ-L 
and experimental resultsO The SU(7V) and 0{N) theo- 
ries are found to be qualitatively correct. In the quantum 
Hall system the magnetic degrees of freedom are carried 
by conduction electrons with S =p-|l/2. Thus single-ion 
anisotropies do not play any role.o 

We describe the anisotropic 2D ferromagnet by a 
Heisenberg model on a 2D square lattice, augmented by 
second- and fourth-order single-ion anisotropies, 
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(ij) 
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where S • S = S{S + 1) and the sum over {ij) is 
over nearest neighbors, counting each bond once. J is 
the exchange constant and B is the external magnetic 
field. Positive anisotropy parameters K2,K^ > fa- 
vor a perpendicular magnetization (easy-axis magnet), 
while the opposite sign favors an in-plane magnetiza- 
tion (easy-plane magnet). Two important special cases 

are = 0, corresponding to uniaxial symmetry, and 
K2 = 0, kI = K^, corresponding to cubic symmetry. 
Equation (H) is viewed as a microscopic Hamiltonian for 
a ferromagnetic monolayer, where the anisotropy param- 
eters are obtained, for example, from first-principles cal- 
culations. The same Hamiltonian can also be used ap- 
proximately as an effective model for films with several 
layers, in which case the anisotropy constants are to be 
understood as averaged parameters containing both bulk 
and surface contributions. Since anisotropy energies (and 
also the exchange coupling J) can be quite different for 
the surface layers and for the bulk, the effective parame- 
ters can be tuned by changing the film thickness. 

Equation (|^) does not include the magnetic dipolar 
interaction 



Si ■ Sj 3 (Si • rij){Sj ■ Tij) 



, (6) 



where = — rj is the separation vector of sites i and 
j, uj — /iog^/iB/47r, /io is the permeability of free space, 
g is the Lande factor, and /^s is the Bohr magneton. 
The dipolar interaction is difficult to treat in the present 
framework, mainly because of its long-range nature. One 
may incorporate it in a simple mean-field way by replac- 
ing one of the two spin operators in Eq. (^ by its thermal 
average, which is, in our notation, the magnetization M. 
In a continuum approximation this leads to 

Hdip = ^ {M,,My, -2M,) • ^ S„ (7) 

i 

where a is the lattice constant. Equation (Q) describes 
the interaction with an effective demagnetizing field. 
This term explicitly breaks spin symmetry in that it 
prefers the magnetization to lie in the xy plane, as is well 
known. In this approximation, the rotational symmetry 
within the plane is retained. A better description would 
also break this residual symmetry. This is the reason why 
a magnetic monolayer has a finite in-plane magnetiza- 
tion even if only the exchange coupling and the magnetic 
dipole interaction are consideredqaS] has been shown us- 
ing a Green's-function approach.BEZI Since we are mainly 
interested in the effects of the single- ion anisotropy terms 
in the Hamiltonian and in spontaneous symmetry break- 
ing, we neglect H^ip for most of this paper. 

We now give an overview over the derivation: After 
replacing the spin operators in the Hamiltonian (|^) by 
bosons, the bosonic system is generalized from two to 
N boson species. We mention that we take care to pre- 
serve the quantum properties of the spin operators in the 
anisotropy terms as we perform the limit — > 00. Then 
the partition function Z is written as a functional integral 
over complex auxilliary fields. The constraint on boson 
number is incorporated using a Lagrange-multiplier field 
and the exchange interaction as well as the anisotropy 
terms are decoupled by means of a Hubbard-Stratonovich 
transformation. Within the mean-field approximation 
the auxilliary fields are then replaced by constants, which 
have to be obtained from saddle-point equations. 

The resulting mean-field free energy Fq depends on the 
four auxilliary parameters P and A. One has to find the 
saddle points of Fq with respect to these parameters and 
check their stability, as discussed in App. ^ A can be 
interpreted as a chemical potential that globally enforces 
the constraint on boson number, whereas the vector P 
is related to the magnetizatioii_see below. Introducing 
rescaled anisotropy parametersEj 
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and the inverse temperature (5 = l/fc^T, the free energy 
reads 



2 ^ 2 V^^^v 
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+ dilnfl-e-^-^^'/2 



(10) 



with the dilogarithm function diln and the effective mag- 
netic field 
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(12) 

(13) 



We emphasize that the whole mathematical apparatus of 
SU(iV) SBMF theory is contained in Eqs. (|l^) -(|l|)pNote 
that these expressions reduce to the isotropic caseQu for 
K2 = = k\ = and that the second-order (fourth- 
order) anisotropics only enter the results if S" > 1 (5 > 2). 
The mean-field magnetization is given by 
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where B' — \B'\. As shown in App. ^ P,; as obtained 
from the saddle-point equations equals the magnetization 
component Mo,i, if the free energy Fq depends on Pi. On 
the other hand, if Pq does not depend on P^, the latter 
quantity can assume any value, while the magnetization 
Mo,i still has a well-defined value given by Eq. (|lj). 

Within the simplified treatment of the dipolar interac- 
tion discussed above, a demagnetizing field proportional 
to (Mo,:c,Mo,y,-2Mo,z) would be added to B'. By us- 
ing a Green's-function method it can be shown that this 
yields an acceptable result for the magnetization if the 
magnetic dipole coupling is small as compared to the ex- 
change coupling.E3 

Finally, we give the results of 0(-/V) SBMF theory with- 
out derivation. The free energy per boson reads 
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p; = p, + i08ii:|p', 

P' =P,„ + 108if|p' 



P; =P, + 6i^2P. 

and the magnetization is 
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(16) 
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with the effective field 



The general form of these expressions is similar to the 
SU(A'^) case and we thus expect qualitatively similar re- 
sults. 



III. RESULTS AND DISCUSSION 

In this section we discuss a number of results obtained 
from SBMF theory and compare them with other ap- 
proaches. Our aim is mainly to show that SBMF the- 
ory yields qualitatively correct results in the presence 
of single-ion anisotropics. We further show that these 
results are comparable to Green's-function methods, al- 
though they are much easier to obtain. We first discuss 
the magnetization Mq and then effective anisotropics. 



A. Magnetization 

SBMF theory correctly predicts the magnetization of 
an isotropic 2D Heisenberg ferromagnet to vanish at fi- 
nite temperatures in the absence of an external mag- 
netic field.Q However, an easy-axis uniaxial anisotropy 
K2 > explicitly breaks the SU(2) spin symmetry down 
to a residual Ising symmetry. A discrete symmetry 
can be spontaneously broken and one expects a finite- 
temperature phase transition between an ordered and a 
disordered phase. Both SU(Af) and 0{N) theory indeed 
yield a finite Tc- Figure |l| shows the magnetization of 
a spin-1 system with a second-order uniaxial anisotropy 
K2/ J — 0.01 in the absence of an external magnetic field, 
as a function of the reduced temperature T /Tc. As usual, 
simple mean-field theory overestimates the magnetiza- 
tion. In particular, at low temperatures the magneti- 
zation approaches M(r — 0) exponentially due to the 
neglect of spin waves. On the other hand, SU(7V) and 
0(iV) Schwinger boson theories and many-body Green's- 
function methodg22l correctly predict an almost linear be- 
havior. However, the absolute temperature scales are 
rather different, as shown in the inset of Fig. |^. Numeri- 
cal results, e.g., from Quantum-Monte-Carlo simulations, 
would be very useful to check the absolute value of Tc. 

One can obtain analytical expressions for by ex- 
panding the saddle-point equations for a small magneti- 
zation Mq.z. The results ard£3 
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for SU(A^) theory and 

rpO{N) _ 



3\ii[l + TTJS^/K2{S -1/2)] 



(20) 



(21) 



for O(A^) theory. Both expressions have the same func- 
tional form. In particular, Tc vanishes in the limit 
K2 0+, as it should, and = for 5 = 1/2. For 
comparison, using a continuum approximation similar to 
the one employed here, the many-body Green's-function 
method yields 
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RPA 



4ttJS{S+1) 



3 ln[l + 37r2 75/4X2(5 - 1/2)] 



(22) 



within the Anderson-Callen decoupling of the second- 
order anisotropyEll Note that all three expressions show 
the logarithmic dependence of Tc on K2, while the de- 
pendence on the spin quantum number differs. 

Since the results of SU(-/V) and 0(iV) theory qualita- 
tively agree with each other, we concentrate on SU(Af) 
in the following. In the case of a hard axis, K2 < 0, and 
kI = 0, the residual symmetry is of XY type. Because 
of the continuous dcscncracy there should be no sponta- 
neous magnctization.l3 This is indeed found within SBMF 
theory. If we include a dipolar interaction by means of 
Eq. (R), a finite in-plane magnetization emerges, since in 
this approximation the dipole coupling acts similarly to 
an external field. 

In the presence of an external magnetic field B one 
always expects a finite magnetization. For an isotropic 
film this magnetization is parallel to the applied field, 
whereas in the presence of anisotropy the magnetiza- 
tion may be oriented along a different direction. In 
addition, a magnetic reorientation transition can occur, 
which may be controlled by varying the temperature T, 
the external magnetic field B, or the film thickness. In 
Fig. |(a) the magnetization components parallel and per- 
pendicular to the magnetic field are shown as functions 
of its strength. The field is applied perpendicularly to 
the easy axis, which is realized by a second-order uni- 
axial anisotropy K2/J = 0.01. Within our theory we 
find a linear increase of the in-plane magnetization along 
the magnetic-field direction below a critical reorientation 
field -Brco- Since the modulus of the magnetization is 
almost constant, the perpendicular magnetization com- 
ponent decreases correspondingly. At B — Brco the lat- 
ter component vanishes continuously, corresponding to a 
second-order phase transition. For B > B^co the magne- 
tization is parallel to the field and increases only weakly 
with B. Figure ^(b) shows a similar situation, with the 
easy axis realized by the fourth-order uniaxial anisotropy 
> 0, while K2 = Kf ~ 0. Again a reorientation tran- 
sition is obtained as a function of the in-plane magnetic 
field. However, since here the magnetization changes dis- 
continuously, the transition is of first order. We expect 



that the dependence of the order of the transition on the 
order of the uniaxial anisotropy is a universal feature, 
since the sanne.. behavior has been observed within sim- 
pler theories. E^la 

To show the interplay of the various interaction terms, 
we now discuss several phase diagrams. For compar- 
ison. Fig. 1^ shows the well-known phase diagram in 



the JC2-IC4; plane for fCl = and B = derived un- 



der the assumption that the magnetization hnhavq s like 
a macroscopic, classical magnetic moment .c§ll3cj The 
anisotropics should be viewed as effective quantities, 
which depend on temperature etc. Roughly, a perpeiu 
dicular magnetization is found for /C2 > and ]Cj- > 0,E£| 
an in-plane one for /C2 < O and ICj; < O, and a canted 
one for /C2 > O and /Cf < O. These magnetic phases are 
separated by continuous (second order) or discontinuous 
(first order) phase transition boundaries, which merge 
in a tricritical point. A coexistence region between the 
perpendicular and the in-plane magnetization is located 
close to the first-order phase transition line. As usual, the 
first-order phase boundary has been drawn where the free 
energies of both phases are equal. 

The simple phase diagram of Fig. ^ does not apply to 
2D systems, since there should be no net in-plane mag- 
netization if solely uniaxial anisotropics are taken into 
account. The reason is that in this case a ground state 
with finite in-plane magnetization would be continuously 
degenerated The origin of this problem of course lies in 
the neglect of fluctuations. 

This picture is changed if fluctuations are included by 
means of SU(7V) SBMF theory: Figure | shows the phase 
diagram in the K2-Kj^ plane for a thic film at a finite 
temperature T/J = 'i for spin S — 2j£3 The anisotropy 
parameters are now the microscopic ones contained in the 
Hamiltonian (||). Note that the scale of the anisotropics 
in Fig. ^ is small, since the temperature T/ J — iis signif- 
icantly lower than the typical temperature scale Ai:JS^ 
of Eq. (|2^) . This means that very small anisotropics are 
sufficient to induce a finite magnetization, since K2 en- 
ters logarithmically in Tc, Eq. (20). 

In Fig. ^ only two phases are present: A phase with- 
out net magnetization located approximately in the re- 
gion K2 < and < and a phase with perpendicular 
magnetization otherwise. These two phases are separated 
by either a first-order or a second-order transition. The 
tricritical point, at which the two merge, has shifted to 
finite, temperature-dependend values of the anisotropics 
K2C > and K^^ > 0. For smaller anisotropics thermal 
fiuctuations overcome the tendency to order. A coex- 
istence region is found in the vicinity of the first-order 
transition. The first-order line is no longer straight and 
approaches the tricritical point with a vertical tangent. 
For negative the second-order transition line is given 
by K2 = const, similar to Fig. ||. Thus the location of 
this transition depends only on K2. At this transition 
the zero-magnetization (M = 0) saddle point becomes 
unstable. The free energy for small M, and thus the sta- 
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bility of this saddle point, is solely determined by the 
lowest-order anisotropy K2. 

In-plane and canted magnetic phases do not appear 
in Fig. ^ since this system is either in the Ising or 
in the XY universality class. However, such phases 
do exist if the in-plane symmetry is explicitely broken, 
e.g., by a single- ion in- plane anisotropy k\ 7^ 0, by 
an external magnetic field, or by the dipole interaction. 
Here we consider a small fourth-order in-plane anisotropy 
k\/J « 4.3 X 10~^ > 0. The other parameters are un- 
changed. The resulting phase diagram is shown in Fig. 
exhibiting now three phases similar to Fig. ^. How- 
ever, the transition between the canted and perpendicu- 
lar phases is now of first order (the first-order line extends 
to arbitrarily large negative K^). The first-order char- 
acter is typically weak — the jump in the magnetization 
is rather small. The coexistence regions are also shown 
in Fig. ^. Note that for even smaller k\ thermal fluctua- 
tions would destroy the in-plane magnetization, resulting 
in a phase diagram similar to Fig. ]^ 

Considering now a thin film with a finite number n of 
atomic layers, the anisotropies have to be viewed as av- 
erages over all layers and can be tuned by varying n. A 
variation of the film thickness thus refers to a certain tra- 
jectory ("anisotropy flow"cj) in the corresponding phase 
diagram, e.g., Fig. |[ In this way we can in principle 
describe reorientation transitions as a function of film 
thickness. 

Finally, we consider the effect of a finite magnetic 
field. We keep the absolute value of the field constant, 
B/ J = 0.01, and change only its direction. For simplic- 
ity we restrict ourselves to the case k\ = 0, and consider 
the case 5 = 2 and T/ J — 1 as an example. The mag- 
netic field is rotated in the xz plane from 9 — to 9 = it, 
where 9 is the angle between the field and the film nor- 
mal. Figure ^(a) shows the phase diagram for the case 
— 0. The magnetization tries to be aligned parallel 
to the magnetic field, but is hindered by the anisotropies. 
For a sufficiently small K2 the magnetization follows the 
field continuously, whereas for K2 larger than a critical 
value it exhibits a jump. This behavior is illustrated 
by the inset in Fig. ^(a), which shows the perpendicular 
magnetization component as a function of 9 for K2 = 0, 
K2 = K2C, and K2 = '2K2c. A negative value of Re- 
does not change this behavior, except that the coexis- 
tence regions becomes more narrow. For Ki- > 0, how- 
ever, the first-order line extends to smaller K2 and splits 
into a two- prong fork, as shown in Fig. |^(b) . Each prong 
ends in a critical point. The metastability regions are 
also indicated in Fig. ^(b). The fork-like feature grows 
for increasing and the critical points even enter the 
negative K2 region for sufficiently large K^. As a side 
remark, in the language of catastrophy theory the phase 
diagram in the {K2, ,9) space shown in Fig. |6| can 
be described by a butterfly catastcpphy. The special cut 
Kj: = 0, Fig. ||(a), shows a cuspEn 



B. Effective anisotropies 

The effective magnetic anisotropies are most generally 
defined through the dependence of the free energy on the 
magnetization direction. The anisotropic part of the free 
energy F{T, 9, (j>) [9 and being the polar and azimutal 
angles of the magnetization) depends on the symmetry of 
the lattice and is usually written with the help of spher- 
ical harmonics YijS:,,^^ or as a series in powers of the 
direction cosines. e3 For example, we give the expression 
for tetragonal symmetry, including also the Zeeman term, 

F = -B • M - /C2 cos^ 61 - /C3 cos^ 9 ~ ICi cos'* 9 

-/c|sin'*0(3/4-|-cos4</>)-... (23) 

The effective anisotropies IC2{T) etc. depend strongly on 
temperature. They have to be distinguished from the mi- 
croscopic anisotropies K2 etc. appearing in the Hamil- 
tonian (||) . Experimentally, the effective anisotropies are 
determined from, e.g., ferromagnetic resonance (FMR) 
experiments. To compare the microscopic anisotropies 
calculated at T = with measurements of the effective 
anisotropies performed at finite temperatures, the tem- 
perature dependence of the latter has to be known. This 
task will be pursued in the present section by use of the 
Schwinger boson technique. 

The main source of the temperature dependence of 
the effective anisotropies is the decreasing magnetization 
M(T).l3 Other sources, such as the population change p£ 
spin-orbit splitted energy levels near the Fermi energy,Lj 
will not be considered here, since the thermal variation 
due to these effects occurs mainly at higher tempera- 
tures. Note that the effective anisotropies depiead on mi- 
croscopic anisotropy terms of various orders ,CJ and can 
even be present at finite temperatures if correspopding 
microscopic ones are absent from the Hamiltoniancj 

For the determination of the effective anisotropies we 
proceed similarly as in Rcf. |42| In order to turn the mag- 
netization into a direction (9, (j)) different from the equi- 
librium one, we apply an auxiliary magnetic field B. For 
the magnetization we use the SBMF result Mq, Eq. (|l4|). 
There are infinitely many choices of the magnetic field 
B that result in the same magnetization direction. This 
freedom of choice introduces some arbitrariness, since the 
resulting effective anisotropies depend on the choice of B. 
We here choose the field B-|= (i?,0,0) along the x axis 
and vary only its strength.^ We insert the SBMF result 
for the free energy Fq, Eq. (|l0|), into Eq. (|2|). To obtain 
the effective anisotropies we calculate the angles 9, (j> and 
the free energy parametrically for a sufficient number of 
values of the magnetic-field strength B and fit Eq. ( |2^ ) 
to the resulting data. 

We now apply this method to a system with uniaxial 
easy-axis anisotropy. In this case the free energy does not 
depend on the angle in the absence of a magnetic field 
and hence the effective anisotropies involving </) vanish. 
Then the free energy has the form 
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i^o + B • M = - 2^ /C„ cos" I 

n=2 



(24) 



We calculate the free energy Fq for equally spaced values 
of the direction cosine cos 6 by varying the field strength 
B. If a first-order reorientation transition takes place 
at some field value B^-^o we use only fields B < B^^o- 
We finally obtain the effective anisotropies by a least- 
square fit. Figure shows the effective anisotropies 
for a thin film with S = 2, J = 100, and uniaxial 
anisotropics K2/J = 0.02, K^/J = 0.01 (filled sym- 
bols) and K2/J = 0.02, K;t/J = -0-01 (open sym- 
bols). We here consider relatively large anisotropies for 
illustrative purposes so that higher-order anisotropies do 
not vanish in the numerical errors. The odd effective 
anisotropies IC3 etc. vanish because of symmetry. A small 
sixth-order effective anisotropy ICq{T) appears at finite 
temperatures, as well as anisotropies of even higher or- 
der, which are much smaller, although no corresponding 
terms are present in the Hamiltonian. The exact val- 
ues of /C2,4(0) for T ^ are /C2(0) = K2 3(3 - 1/2) = 
K23^ aiidJfA{0) = Ki S{S - 1/2){S - l){3 - 3/2) = 
^±5.4 ggg (|) (|) Q^j. approach re- 

produces these limiting values. Note that the special 
generalization of the anisotropy terms from SU(2) to 
SU(iV) discussed in App. ^ is required to obtain these 
results. For T Tq-^he effective anisotropies vanish 
as /C„(T) cx M"(T)MeI Therefore, lCi{T) decreases 
faster with increasing temperature than /C2(T) and ex- 
hibits a more smooth behavior near Tc- The effective 
anisotropies shown in Fig. ^ qualiiatively agree with 
many-body Green's-function results.E3 

Finally, we emphasize that the calculation of effective 
anisotropies with this method suffers from ambiguities 
due to the freedom of choice of the magnetic field and 
the particular fitting procedure employed. Furthermore, 
this derivation employs a static calculation of the mag- 
netization and free energy in order to describe a dynam- 
ical (FMR) experiment. These problems are common to 
all calculations of effective anisotropies we are aware of. 
A full quantum-mechanical, dynamical theory of FMR 
would be very useful. One possible approach would be 
to incorporate spin dynamics into the Schwinger boson 
theory. 



IV. CONCLUSIONS AND OUTLOOK 

In summary, we have presented the Schwinger bo- 
son mean-field theory for anisotropic 2D ferromagnets, 
using two different bosonizations. In this study we 
have focused on the theoretical treatment of single-ion 
anisotropies of a Heisenberg-type Hamiltonian. The 
anisotropy energies break the full SU(2) spin symme- 
try, leading to a lower residual symmetry, e.g., of Ising 
or XY type. We have shown that the results for the 
magnetization satisfy the Mermin- Wagner theoremB also 



in these cases, in that a finite spontaneous magnetiza- 
tion is only found if the ground state is not continu- 
ously degenerate. The theory can describe the various 
magnetic order-disorder transitions as well as the reori- 
entation transitions found in experiment. The bosoniza- 
tion method can serve as an alternative technique for 
quantum spin systems, besides the many-body Green's- 
function approachHO Our results for magnetization and 
effective anisotropies qualitatively agree with the ones 
obtained by the latter method. Microscopic calcula- 
tions of effective anisotropies as functions of tempera- 
ture are important, since they — and not the microscopic 
anisotropies appearing in the Hamiltonian — are obtained 
from experiments. However, the absolute value of the 
critical temperature at which the spontaneous magneti- 
zation becomes finite strongly depends on the details of 
the theory. Compared to the Green's-function method, 
Schwinger boson theory has the advantage of treating 
any spin value S and any orientation of the magneti- 
zation relative to the anisotropy axes and the external 
magnetic field on the same footing. It is also numeri- 
cally much less demanding. Other types of anisotropies, 
for example of hexagonal surfaces, can be incorporated 
in a straightforward manner. Finally, Schwinger boson 
mean-field theory is a well-controlled approximation in 
the sense that fluctuations about the equilibrium solutioii 
can be treated systematically within a 1/A^ expansion.l3 
We conclude with enumerating possible generalizations 
of the theory that are of particular interest for a realis- 
tic description of magnetic thin films: The first is the 
explicit description of thin films of several layers, where 
the exchange couplings, the magnetic moments, and the 
anisotropies are microscopic quantities which depend on 
the layer index. This can be done by a straightfor- 
ward generalization of the theory, where one only has 
to keep track of additional layer indices. The second 
is the inclusion of the long-range magnetic dipole in- 
teraction. This interaction is an additional source of 
anisotropy, which acts quite differently from the single- 
ion anisotropies, resulting for instance in the formation 
of magnetic domains O We reiterate that both the ef- 
fect of several layers and the dipole interaction can be 
treated in the present theory as it is, albeit only in a 
mean-field manner. The third interesting generalization 
is the inclusion of magnetic dynamics. This would allow 
to describe the precession of a non-equilibrium magneti- 
zation and constitute a microscopic theory of ferromag- 
netic resonance. We leave these questions as projects for 
the future. 
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APPENDIX A: 

In this appendix we present the Schwinger boson the- 
ory for anisotropic ferromagnetic films in some more de- 
tail. To derive the mean-field equations we start from 
the Hamiltonian (|^). We replace the spin operators at 
each lattice site by two Bose fields 6|i and 6|i according 
to Eqs. The fields b^i and bii have to satisfy the 

constraint at each lattice site. With the matrix 



S(i) 



(Al) 



the Hamiltonian reads, up to an additive constant, 



H = - 



B 



2 ^ 



^ E + (0]' ) (A2) 



Ae = 



[S^ -{S- 2lNf] = 2K2,nS - 



N 



(A4) 



is kept constant for all N. In particular, Ae = for 
5 = 1/2 and all N. This is achieved by setting 



J^2,N = — 7777^2- 



(A5) 



S-l/N 

In the mean-field approximation, N ^ oo, we get@ 

ih = K2,oo ^K2{S~ 1/2)/ S. (A6) 

Analogously, the A^-dependent fourth-order anisotropics 
are obtained. 



(S - 1/N){S - 2/N){S - 3/A) 



K, 



4 = Kt:^ = Ki-^^ {S l/2)(5 -1){S- i/2)/S\ 

(A8) 



with the Pauli matrices cr = (a^, tr^, ct^). Summation 
over repeated greek indices is implied. 

This expression is now generalized to N Bose fields ba ■ 
Then S becomes an A x A matrix of operators b]j)p and 
the constraint is generalized to 



bib^ = NS. 



(A3) 



The Pauli matrices are generalized to N x N block- 
diagonal matrices with the original Pauli matrices re- 
peated along the diagonal. We keep the notation etc. 

There is some freedom in the generalization of the 
Hamiltonian to N boson flavors — we mainly have to make 
sure that the correct N = 2 case is recovered. If we 
naively generalize the second-order anisotropy term as 
-{K2/2N) (0-^)^3 'S'f]^, we run into problems, which 
are most easily illustrated for the case S = 1/2. For 
A = 2 the operator has the eigenvalues ±1/2 so that 
(S"^)^ only has the eigenvalue 1/4. This is, of course, the 
reason why the K2 term is irrelevant for S = 1/2. How- 
ever, for S" = 1/2 and A > 2 (A even) we find A/2 + 1 
different eigenvalues of the z spin component {<J^)'pS^ 
and A/4 -I- 1 (or A/4 -1-1/2) different eigenvalues of the 
K2 term for A (not) an integer multiple of 4. Thus in 
the limit A — > 00 the second-order anisotropy term has 
infinitely many different eigenvalues even for S = 1/2, 
which looks like a semi-classical approximation. Clearly, 
the mean- field results would depend on K2 for 5 = 1/2, 
which is incorrect. To remedy this problem we introduce 
an A-dependent anisotropy K2,n such that the energy 
difference between the two lowest energy states of the 
second-order anisotropy term. 



The fourth-order anisotropy contributions now correctly 
vanish for S* < 2 even at the mean-field level. 

In order to write down a functional integral over the 
Bose fields, we have to normal order the bosonic oper- 
ators. The exchange term is trivial, since the operators 
at different lattices sites commute. By commuting the 
operators in the anisotropy terms, using the constraint 
^), and dropping a constant we obtain 

HsuiN) = -^E^^"«^"W - ?E'^;^^"« 



K' 



2.N 



2N 

2A3 



i 



K 



2^E^([(-^)?^^«]' + [(^')^^^w]' 



(A9) 



where : : denote normal ordering. The second-order 
anisotropy has obtained a contribution from commuting 
the bosons . in , t he fourth-order terms, similar to other 
approaches ,l3l3 



K2.N = K2.N + 



6A5 + 4 
A2 



6NS + 4 ^\ 



A2 



4,N- 



(AlO) 



If we want to do a consistent theory to order zero in 
1/A we can just drop the corrections. However, since 
the prefactors {6NS + 4) /A^ = 35" -I- 1 are not small for 
A = 2 it may be better to keep them. Note also that 
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they cancel for cubic symmetry. In the following we omit 
the prime at K2 but keep these points in mi*id. 

Next, we go over to the continuum limitEl and write 
down the functi onal integral for the partition function. 
The constraint (A3) is implemented using a Lagrange 
multiplier A(r, r) at each point in space r and imaginary 
time T.Q A is to be integrated from —ioo to +ioo but we 
can deform the path of integration between these lim- 
its. The path of integration has to cross the real axis on 
the positive side to avoid a branch cut. The partition 
function reads 

Z = j D^ba,DXexp(^-^J^ '^dT J (frC[b;X]j , (All) 

where now ^^(r, r) is a complex field. The Lagrangian is 
C[b; A] = ^b^drb^ + JSid,bl)id,b^) 
L bl{d,b})bp{d,b^) ~ ^CT^b^b^ 



N 
K2N 



2Na? 



2a2 



K. 



\b*J)a 



NSX. 



(A12) 



Here, a is the lattice spacing, dr is the time derivative, 
and dj is a spatial derivative along the j direction. Sum- 
mation over repeated indices is implied. 

In the next step, the terms containing more than 
two ba fields are decoupled using Hubbard-Stratonovich 
transformations, thereby introducing additional auxiliary 
fields Q, P, and R, which have to be integrated over in 
the functional integral. In a short-hand notatign the de- 
coupling of the exchange term may be writterO 

b:{d,b*p)bf^{d,b^) ^ NJQ.Q, + iJQ, bl{d,b^) 

- iJQj {djbl)b^. (A13) 

The K;^ term contains eight Bose operators, which are 
decoupled into four according to 



K, 



NK 



N 



(A14) 



To make the transformation well defined, Rz is to be 
integrated from —00 to +00 if > but from —ioo 
to -|-ioo if < 0. The remaining four-boson term is 
grouped together with the K2 term and decoupled. 



2,N 



2N 

N{K2M 



{K2,N + 2Klj^Rz)Pz{a' (A15) 



Here, the direction of integration of must be chosen 
such that [K2.N + "^K^j^RzjPz is real and positive at 

the end points. The decoupling of the kI terms is done 
similarly, introducing additional fields Rx,y and Px,y 

Up to this point no approximation, apart from the con- 
tinuum limit, has been made. Now we consider the mean- 
field approximation for this system. This is exact for 
N —^ 00 but serves as a lowest-order approximation for 
N = 2. The mean- field theory is defined by all auxiliary 
fields being constant in time. Since the system is ferro- 
magnetic, we assume that they are constant in space as 
well and write them as A, Q, P, and R. Then, the action 
is bilinear in the fields Written in terms of Fourier 
transforms the mean-field partition function reads 

Zo = Z, j Yl d%^{k,iLu„) expf- y" d^kY,^o[b]] 

(A16) 

with the 6ct-independent part 

= exp I - MNpjq, ■ Qa2 + Afa^NSfiX 



-MN^-^Rl 



-MN^-^Pl 



R.r. + R, 



-MN(3Kl X 



RxP. 



RyPy 



(A17) 



where M is the number of sites. The Lagrangian is 



-l-a Aj b*^{^^,^Un)ba{^^,^uJn) 
+ 2kiRzPz{a'')1 

+ 2kl [RxPx{o^rp+RyPv{oy)f\^ b}ba. (A18) 

After integrating out the ba fields and performing the 
sum over field indices a, we obtain, up to an irrelevant 
constant factor. 



Zc exp^ - y ^^^y d^fc^lndetM(k, 



lU!,-, 



(AI9) 



with the 2x2 matrix 



M = {-ifihujr, + (iJSk^a^ + A)l - • cr 

- (ik2 Pza'-~2 13 RzPza' 

- 2 pkl (RxP^a'' + RyPyay). 



(A20) 
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Note that Zq only depends on A and Q through 



i?,; = P. 



(A29) 



A = a2/3A - ^ Q • Qa^ 



(A21) 



just as in the isotropic case.0 Evaluat ion o f the frequency 
sum and momentum integral in Eq. ( A19 ) yields the free 
energy per boson 



1 



-S'A4 



InZn 



— - — r ■, 



PKi -2 

(-2 , -2N 



( i?,P, + RyPy 



1 



diln 



-h+l3B' /2 



+ diln ( 1 - e-A-/5s72 



with B' = \B'\ and 



By - By 



AKlRyPy, 



B', = B, + 2K2Pz+4KiR,P,. 



(A22) 

(A23) 
(A24) 
(A25) 



The dilogarithm function is defined by the series repre- 
sentation 



dilnx 



E 



(1 -x)" 



(A26) 



We are mainly interested in the magnetization 
„ dFo 1 B' 



In 1 - . 



dB 



-A+f3B'/2\ _ In (x _ ^-A-l3B'/2 



Stt^JS B' 



{A27) 



We have used that the derivatives of Fq with respect to A 
and P vanish at the saddle point. A factor of 2 appears, 
since in the physical system two bosons make up one 
spin. For the isotropic case B' equals B and we recover 
the results of Refs. ^ and 0. 

Finally, we have to determine the auxiliary fields by 
finding the lowest stable free-energy saddle point. We 
set the derivatives of InZp with respect to A, P, and R 
to zero and solve the resulting equations simultaneously. 
From dlnZo/ dA = we obtain the saddle-point equation 



1 



8Trf3JS 
In f 1 - 



In fl - e-A+zas'/a 



^-A-PB'/2 



(A28) 



which can be analytically solved for aE By taking the 
derivatives with respect to the other auxiliary fields one 
finds that for all saddle-point solutions 



This equation either follows directly from the saddle- 
point equations or, for special values of Ri, the free en- 
ergy is found to be independent of the corresponding Pi 
so that we can choose Pi to satisfy Eq. ( A29 ). We end up 
with a free energy which depends on the four auxiliary 
parameters A and P, 



pFo-- 



2 
1 



diln 1 



-A+/3S72 



8ttI3JS 
+ diln { 1 - 6-^-/5^72 



with 



B'^ = B.,,+AKlPl 
B'y = By + 4Klpl, 
B'^B,+2K2P,^ 



4KtP, 



(A30) 

(A31) 
(A32) 
(A33) 



In the uniaxial case Kl = there are only two relevant 
parameters, A and P^- 

By solving the saddle-point equations d In Zq /dPi = 
we find that there are typically two solutions for each 
component Pi. One is Pi — Mo^i, i.e., Pi is the cor- 
responding magnetization component, and the other is 
given by B'^ = Bi. For example, in the uniaxial case, 



kI = {KH = 0) the free energy does not depend on P^ 
and Py whereas Mq^x and Mo^y have well-defined val- 
ues, which need not be zero in the presence of an ex- 
ternal field. Thus Px.y cannot equal Mo^x:,y The saddle- 
point equations are solved by the second type of solution, 
B'xy = Bx^y, since K^- = 0. On the other hand, the z 
component is Pz = Mq.z. 

It turns out that the relevant saddle point of the free 
energy Fq is not a minimum with respect to all four pa- 
rameters. For example, it is always a maximum with re- 
spect to A (or A). Intuitively, one would expect A t o run 
off to infinity in this case. However, since A in Eq. ( All ) 
has to be integrated from — ioo to -|-ioo, fluctuations in A 
are along the imaginary direction, A(r, r) = A-|-iAA(r, t), 
where AA is real. In particular, this holds for global 
changes of A. With respect to AA, the free energy in- 
deed has a minimum at A A — 0. For the Hubbard- 
Stratonovich fields P and R the direction of fluctuations 
depends on the values of the anisotropy constants. The 
stability of each saddle point against global changes of the 
auxiliary parameters has to be checked in order to find 
the lowest stable one. Note that for positive anisotropy 
constants the fluctuations of P and Q are purely real and 
the stability analysis becomes straightforward. 
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FIG. 1. Magnetization as predicted by various theories as a 
function of reduced temperature T/Tc for a film in a vanishing 
external magnetic field, B = 0, with an easy axis defined by 
K2/J ~ 0.01, spin S — 1, and no higher-order anisotropies. 
The magnetic dipole coupling is neglected. The inset shows 
the magnetization as a function of the absolute temperature 
T/J in units of exchange J. The line denoted by 'RPA' 
refers to a calculation within a many-body Green's-function 
approach uith Anderson-Callen decoupling of the single-ion 
anisotropy.Ej 
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FIG. 2. Reorientation transitions as a function of external 
magnetic field for a film with exchange J = 100 and spin 
S = 2 at temperature T/J = 1. (a) Film with second-order 
easy-axis anisotropy K2/J = 0.01 and K^r ~ = 0. The 
field is applied in the x direction, i.e., perpendicularly to the 
easy axis. The solid and dashed lines represent the x and z 
components of the magnetization. In this case a second-order 
reorientation transition takes place, (b) The same quantities 
for a film with forth-order easy-axis anisotropy / J = 0.02 
and K2 = kI = 0. Here, a first-order reorientation transition 
takes place, denoted by the arrows. The dotted lines refer to 
meta-stable states of higher energy. 
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FIG. 3. Phase diagram for the effective anisotropies K,2{T) 
-d /C£(2ZLjEpr a vanishing external magnetic field und 
0.c3'E3o Three different phases are found, defined by 
different directions of the magnetization with respect to the 
film plane (perpendicular, in-plane, canted), indicated by ar- 
rows. These phases are separated either by a first order (solid 
line) or by a second order phase transition (dashed lines), 
which merge at the tricritical point K.2c ~ IC4C ~ 0, corre- 
sponding to the isotropic, SU(2)-symmetric case. A coexis- 
tence region between the perpendicular and in-plane phases 
exists near the first order phase transition boundary, indicated 
by the dotted lines. 
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FIG. 4. Phase diagram of a thin film in the K2-K4 plane in 
a vanishing external magnetic field as calculated from SBMF 
theory. The parameters are spin 5 = 2, exchange J = 100, 
temperature T/J — 3, and in-plane anisotropy ifj = 0. The 
phase to the lower left of the phase diagram now exhibits 
no net magnetization, M = 0. The open circle denotes the 
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FIG. 5. Phase diagram for a thin film as in Fig. ^ but with 
a finite in-plane magnetic anisotropy K\jJ ~ 4.3 x 10~*, 
big enough to allow a finite in-plane magnetization at this 
temperature. The square denotes the point of full cubic sym- 
metry {K2 ~ 0, — K\). The dotted lines are borders of 
two coexistence regions. In the left one both perpendicular 
and in-plane magnetization are metastable, while in the right 
one perpendicular and canted states are metastable. 
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FIG. 6. Phase diagram for a thin film in the K2-O plane 
as calculated from SBMF theory. An easy-axis magnet is as- 
sumed, with K2 > 0. A magnetic field of strength B / J = 0.01 
is applied in the xz plane, forming an angle 6 with the film 
normal. We assume J = 100, S = 2,T = J, and kI = 0. Fur- 
thermore, in (a) Kj- = 0, whereas in (b) /J « 5.3 x 10~^. 
The thick solid lines denote first-order transitions ending in 
critical points. The dotted lines are the boundaries of co- 
existence regions. The inset shows the magnetization as a 
function of 6 for the case (a) with K2 = 0, if 2 = if 2c, and 

K2 = 2A-2C, if 2c. 
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FIG. 7. Effective anisotropies /C2, ICi, and ICe as functions 
of temperature calculated from SBMF theory. A thin film 
is assumed with S = 2, J = 100, K2/J = 0.01, kI = 0, 
and Ki/J = 10"^ (filled symbols) and K^/J = -10"^ 
(open symbols). We use rather large microscopic anisotropies 
to demonstrate more clearly the appearance of a small 
sixth-order effective anisotropy /Ce, which is shown scaled by 
a factor of 50. The big symbols denote the exact values at 
T = 0. 



